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ABSTRACT 

East and effective localization of gravitational wave (GW) events could play a crucial role in identi¬ 
fying possible electromagnetic counterparts, and thereby help usher in an era of GW multi-messenger 
astronomy. We discuss an algorithm for accurate and very low latency (< 1 second) localization of 
GW sources using only the relative times of arrival, relative phases, and relative signal-to-noise ratios 
for pairs of detectors. The algorithm is independent of distances and masses to leading order, and 
can be generalized to all discrete sources detected by ground-based detector networks. Our approach, 
while developed independently, is similar to that of BAYESTAR with a few modifications in the al¬ 
gorithm which result in increased computational efficiency. Eor the LIGO two detector configuration 
(Hanford-I-Livingston) expected in late 2015 we find a median 50% (90%) localization of 143 deg^ 
(558 deg^) for binary neutron stars (for network SNR threshold of 12, corresponding to a horizon dis¬ 
tance of ^ 130 Mpc), consistent with previous findings. We explore the improvement in localization 
resulting from high SNR events, finding that the loudest out of the first 4 (or 10) events reduces the 
median sky localization area by a factor of 1.9 (3.0) for the case of 2 GW detectors, and 2.2 (4.0) 
for 3 detectors. We consider the case of multi-messenger joint detections in both the GW and the 
electromagnetic (EM) spectra. We specifically explore the case of independent, and possibly highly 
uncertain, localizations, showing that the joint localization area is significantly reduced. We also show 
that a prior on the binary inclination, potentially arising from GRB observations, has a negligible ef¬ 
fect on GW localization. Our algorithm is simple, fast, and accurate, and may be of particular utility 
in the development of multi-messenger astronomy. 


1. INTRODUCTION 

One of the most exciting potential scientific payoffs of 
the coming age of gravitational wave detections is multi¬ 
messenger astronomy: the observation of a source in both 
the gravitational wave (GW) and electromagnetic (EM) 
spectrum. Although in some cases the EM sources may 
naturally be observed in coincidence with the GW sig¬ 
nals (e.g., short gamma-ray bursts [GRBs]), the more 
common case is expected to consist of GW detection 
followed by EM exploration of the GW sky localization 


error boxes (Metzger and Berger 2012 Chen and Holz 


20131. The most promising sources for ground-based GW 


instruments are stellar mass binary inspirals and merg¬ 
ers of neutron stars and/or black holes. It is not known 
if these systems are accompanied by EM signatures, al¬ 
though there is growing evidence that short GRBs result 
from binary coalescence ( Koce vski et al. 2010 Eong et al. 
2013 Eong and Berger| 2013(). I'hese systeins may also 
produce N-ray, radio, optical, and/or infrared afterglows 
(e.g, kilonovae p owered by the radioactive decay of r- 
process elements |Berger et al.||2013l Tan^vir et al. 


Freiburghaus et ar T999| Metzger et al. 2010|) ). 


2013 


There 

will be tremendous interest in identifying possible coun¬ 
terparts to LIGO sources. The success of these endeavors 
will depend crucially on the ability of the LIGO network 
to quickly and accuratel y localize GW sources on the sky. 

The advanced LIGO (Harry and LIGO Scientific Col¬ 
laboration 2010) detectors are currently under develop¬ 
ment, and their first science-quality data are expected to 


arrive in Fall 2015 (LIGO Scientific Collaboration et ah, 
2013|) . The envisionM operational schedule consists of an 


initial GW network composed of the two LIGO detectors 
(Hanfo rd [H] and Liv ingston [L]), with the Virgo [V] de¬ 
tector (|Accadia|2012|) joining the network at a later date 
and at lower relative sensitivity. GW source localization 
mainly relies on triangulation based on the arrival times 
of the GWs at the different si tes of the detector net - 
work, as has been discussed in iFairhurstl (120091120111). 


The ability to triangulate a source is crucially dependent 
upon the number, bandwidths, and separations of detec¬ 
tors in the GW network, as well as the signal-to-noise 
ratio (SNR) within each detector. 

There are two basic approaches to localizing a GW 
source. The first, exemplifi ed by the BAYESTAR 


pipeline (Singer and Price 2015) in LALInference, utilizes 
the output of the LIGU detection pipeline to directly es¬ 
timate possible sky localizations (IKasliwal and Nissanke 
2014 Sidery et al.||2014 Singer et al. 2014). I'he detec- 


tion pipeline produces rapid estimates of such quantities 
as time of arrival, phase, SNR, and chirp mass, but with 
fairly broad uncertainties. An alternative method to de¬ 
termine sky localization is to run the f ull parameter esti¬ 
mation pip eline within LALInference (Veitch et al. 

Aasi et al. 2013). This entails fitting for the raw “data 


-|- noise”, utilizing Markov Chain Monte Carlo (MCMC) 
and nested sampling techniques to estimate the relevant 
parameters of the binary. Although this approach gen¬ 
erates the most reliable and tightly constrained values 
of parameters, including the best estimates of the sky 







































































2 


position, it requires significant computational resources 
and can take hours, if not days or longer, to run. The 
performance of BAYESTAR and full parameter estimate 
sampling methods are comparable for GW networks com¬ 
posed of two interferometers, while for three or more de¬ 
tect ore_the_Jull_parameter estimation pipelines are supe¬ 
rior (Singer et al. 2014). This is because, as currently 


configured, the detection pipeline does not produce any 
information about the Virgo data if that detection is sub¬ 
threshold (p < 4). 

In this paper we develop an extremely low latency 
(< 1 second after detection) localization algorithm in¬ 
corporating the timing, phase, and SNR output from the 
detection pipeline to estimate the localization area. We 
take as input the relative timing and phase, and the SNR 
ratios of the GW detectors. From these parameters we 
estimate solely the localization area, marginalizing over 


localization algorithm utilizes the log of the SNR ratios, 
rather than the value of SNR at each detector. This ob¬ 
viates the need to infer the binary chirp mass or distance, 
saving the time (and potential errors) of fitting for these 
parameters, only to marginalize over them afterwards. 
By taking the logarithm, the Gaussianity of the likeli¬ 
hood is approximately preserved. Second, we calculate 
the sky location prior from the antenna power patterns 
and the relative sensitivities of the detectors. This prior 
applies to all discrete sources detected by ground-based 
detector networks, and can be pre-calculated in advance 
of utilization of the pipeline. Finally, our algorithm ex¬ 
plicitly marginalizes over only two unknown parameters: 
the binary orientation and inclination. Since this is a 
small parameter space, we can pre-grid the entire sky and 
pre-calculate this marginalization, allowing our localiza¬ 
tion algorithm to consist primarily of lookup tables and 
summation, rendering it simple and fast. Our alogrithm 
can be run in < 1 second on a laptop, offering the poten¬ 
tial for real-time localization as the GW source evolves 
within the detectors. Since our method has been devel¬ 
oped completely independently of the LIGO code efforts, 
it also provides a consistency check of the existing LIGO 
codebase (and in particular BAYESTAR). 

We use our localization algorithm to explore issues re¬ 
lated to EM follow-up. For example, we consider the 
dependence of the localization area as a function of SNR 
of the events, explicitly determining the expected im¬ 
provement in localization from high-SNR events. The 
incidence of these high-SNR events can be predicted an- 
alytic ally, and depends so lely on the measured event 
rate (Chen and Holz 2014). We show that EM priors 


on the inclination do not significantly help GW localiza¬ 
tion. We also discuss the potential of cross-correlating 
LIGO sources with other all-sky observatories, such as 
Fermi GBM, to dramatically improve sky localization. 

2. LOW LATENCY LOCALIZATION ALGORITHM 

From Bayes’ Theorem, the probability of a source be¬ 
ing located at a sky position {9, ^), given a set of observ¬ 
ables e, can be written as a Bayesian posterior: 




/(e1^,0)/(g,(^) 

/(e1 


( 1 ) 


The prior on source location can be estimated 
from the detector network’s antenna power pat¬ 


tern ( Sathyaprakash and Schut^|2009[ Schutz|[2MI ), 


n{9, (j), L, -Ip) = F^{9, (j), '0)(l-|-cos^ t)^-|-4Fx (9, (j>, 'tjj)cos^ l, 

( 2 ) 

and sensitivity, 


I7 = 


j-7/3 

Mf) 


df, 


( 3 ) 


where ip is the orientation of the binary within the plane 
of sky, L is the inclination angle between the binary’s ro¬ 
tation axis and the line of sight, and Sh{f) is the detec¬ 
tor’s power spectral density. The prior can be calculated 
from: 


all other parameters. Our ap 

proach is very similar to 

BAYESTAR ( 

Singer and Price 

12015 

Singer et al.||2014 

1 where 

with a number of important c 

ifterences: First, our sky 


= / f{9,(p,L,ip)dLd'tp, 


f{9,(p,L,ip)!X. fli (6>, (p, Ip) Ir,}j sin 6» sin (5) 

and where the sum goes over each detector in the network 
{i =H, L, V,...). The likelihood function is expressed 
as a function of the sky location, {9,(p), and the GW 
observables, e: 


( 4 ) 


/(e|6», (()) oc exp (-Ay?(6», 4>)/2) 


( 6 ) 


where Ax| can be calculated from a chi-square minimiza¬ 
tion process: Axl{9,(p) = xl{d,(p) - x|„,in(6'min, <(>min)- 
For our algorithm the observables are the arrival times of 
the signal, U, the phase, rji, and the SNR, pi, measured at 
each detector. If the SNR measurements are independent 
of the arrival times and phases, the chi-square expression 
can be separated: x|(6>, 4>) = X({d, (p) +Xp(^) where C 
represents the contribution from arrival time and phase, 
and p represents the contribution from SNR (see below). 
We assume the errors in arrival time, phase, and SNR 
are Gaussian and independent b etween detectors. This 
appears to be a fair assumption (|Berry et al.||2015|), al¬ 
though our method can be generalized it the error prop¬ 
erties are known. 

The differences in the arrival time between detectors 
allows us to triangulate for the sky position of the source. 
Using only this information, every pair of detectors will 
localize the source to an annulus on the sky. Three detec¬ 
tors with two intersecting annuli will localize the source 
to two intersection regions, while four or more detectors 
can localize the source to a singl e distinct region . This 
technique has been discussed by Fairhurst (2009 2011) 
using a flat-sky approximation; we rederive these results 
in full generality. 

The phase of a binary, p, at time t, measured in the 
ith detector is a function of the detector orientation, 
and the binary sky location, i nclination, and orientation: 
p(t) + rii(9, (p, 6, pj). Following Sathyaprakash and Schutz 
(2009), we find: 


r]i{9,(p,L,ip) = tan 


2F^^i{9,(p,ip)cosi 

F+,j(6», (p, ip){l + cos2 i) 


By considering only the phase difference between detec¬ 
tors, we avoid needing to determine the absolute phase of 
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the binary. The sky location remains degenerate with the 
inclination and orientation, and this can be subsequently 
marginalized over. There are additional phase differences 
due to the separation of the detectors on Earth; this can 
be inferred from the difference in signal arrival times. 
The phase and timing measurements are therefore corre¬ 
lated. 

We follow th e app roach of Grover et al. (2014); 


Fairhurst (2009, 20111) to estimate the timing and phase 


covariance matrlSd These errors depend on the SNR 
and the detector effective bandwidth of the source, cr^ = 

P - 


cov(t, T]) = 


2 2 
2 2 



r 1 7 1 


{27Tpcrf)^ 2-Kp'^cr^ 


f p 

-* 

27I-p2<T2 {p(Tf)^ 


(7) 


/ 1 1-, ( f\ 12 

/"■ df. The error in measured ar- 

^h[f) 

rival times and phases are expected to be distributed as a 
2D Gaussian centered at the true value with covariance 
given by Eq. In reality, the simulated bin ary neu¬ 


tron stars merger signals in Singer et al. (2014) follow a 


Gaussian distribution with a standard deviation that is 
1.4 and 1.3 times larger in arrival time and phase, respec¬ 
tively. The difference in arrival times and phase between 
pairs of detectors is a linear function of the independent 
observables, ti and iji'. 


^tij — ii tj : 

Afjij = f]i - r]j, = H, L,V,... with i^j. 

Thus the measured timing and phase differences can still 
be described by 2D Gaussians. If there are N detectors 
in a given GW network, there are — 1 independent 
timing and phase differences. Without loss of generality 
we can denote one of the detectors as reference detector I, 
and define all time-of-arrival and phase differences with 
respect to this fiducial detector: 

Ati — tj tj, 

Arji = r]i -T]i, I. (8) 

The chi-square value is: 


where 




AU - At^{9,ip) 

Arii - Ar]i{9,4>,L,'ijj) 


and Vij is the covariance matrix. Since all the differ¬ 
ences are in reference to the same detector /, any two 
differences are not independent. 

Following Eq. we find that when i = j the compo¬ 
nents of Vij are given by 


If * 7^ J, 


2 2 

^AtiAtj — ) 

2 2 

^ArjiArjj — 

^AtiArjj — ^AriiAtj (10) 


We marginalize over using Eqs.[^and[^ 

= -2 log y bV') didii^ 

In addition to the timing and phase differences, the rel¬ 
ative values of SNR received by each detector in the GW 
network also provides constraints on the source location. 
Different detectors are sensitive to different polarizations 
of the incoming waves, and the relative response in each 
detector depends upon the binary sky location (0, (j)) and 
orientation well as the detector placement on 

the Earth, orientation of the detector arms, and the de¬ 
tector sensitivity, Sh{f)- The SNR measured in a given 
GW detector depends on the antenna response of the de¬ 
tector, as well as the binary source parameters, such as 
the luminosity distance, redshift, and chirp mass. Aside 
from chirp mass, these binary parameters are poorly con¬ 
strained by the detection pipeline. However, since for 
our purposes we are only concerned with sky position, 
we can eliminate the other binary parameters by taking 
the ratio of the SNR measured at pairs of detectors. This 
ratio is solely a function of the detectors’ antenna power 
patterns (Eq. and sensitivities (Eq. : 


^3 




We estimate the Gaussian uncert ainty in SNR follow¬ 


ing th e Fisher matrix calculation (Cutler and Flanagan 


1994): 


1 


( 12 ) 


which corresponds to a ^ 10% error in SNR measurement 
for a detection threshold of pth ~ 12. The measured SNR 
will be Gaussian distributed about the theoretical value 
p = G(/9, CTp). Unlike the arrival time difference, the SNR 
ratio is not linear in SNR, so the error in the ratio is 
not described by a Gaussian. However, by taking the 
logarithm of the SNR ratio: 


= log — = logp* - iogpj, (13) 
Pi 


we find that the error in R approximates a Gaussian with 
variance 0 



^2 _ ..,2 , 2 

^AtiAti — ^tiU w 

2 — ^2 I 2 

^A7iiAr)i ^TiiTji ^7]iriii 

2 2 2 2 
^AtiArii = ^AriiAti = O’tirji + ^Urn' 

^ If variables are correlated, the variance of their sum is the sum 
of their covariances. 


Similar to the case with arrival times, there are — 1 
independent SNR ratios from a network consisting of 
N independent detectors. Therefore we define Ri = 

^ The error becomes increasingly non-Gaussian when the SNR 
in one detector falls below ~ 5. This leads to a loss of the catch 
percentage at high confidence level, though this deviation is in¬ 
significant (< 1 %). 
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Fig. 1.— Cumulative 50% con fidence localization areas of the 
630 simulated sources taken from |Singer et aL| (|2014|| detected by 
LIGO Hanford+Livingston with the noise curve m TTgo-tl200307 



Fig. 2.— Cumulative 50% co nfidence localization areas of the 
81 simulated sources taken from |Singer et al.|(|2014ll detected by 
LIGO Hanfor d+Livingston and Virgo with tire noise curv es in ligo- 
tl200307 andjLIGO Scientific Collaboration et al.|(|2013ll. We only 
consider events witn biMrt greater tnan 4 m all tnree detectors. 


log(p//pi), i ^ I with an arbitrary reference detector 
/. The chi-square for SNR is then: 


Li 

04) 

where Vij = cov{Ri, Rj) is the covariance matrix. Hi = j 
then 

2 _ 1 1 


while ii i ^ j then we find: 

= \ \^R, + - var(i?i - Rj)] 

_ 1 / 1 1 1 1 1 1 \ ( 15 ) 

~ 2\7i^ 7.^ 7,~ 7.~ 7,) 

1 



Fig. 3. — We compare the improvement of using time of 
arrival+SNR and time o f arrival+SNR+pha se. The simulated 
sources were taken from |Singer et al.| ||2014|) and localized with 
LIGO Hanford and Livingston. 


By taking the logarithm of the power ratio (Eq. 131 we 
are able to estimat e th e covariance matrix analytically. 

As shown in Eq. |14| the directional information is de¬ 
generate with binary inclination and orientation. We 
marginalize over (t,'0) using Eqs. |^and[^ 


Xp(6',0) = -2 log 




(16) 


The sky prior (Eq. also contributes to the poste¬ 
rior, so we define the total chi-square as: X^(6*, </>) = 
x|(0, ^) — 21og/(0,(/>). We grid the sky uniformly in 
(cos 0, (j)) (or in {0, </>) when close to the poles) calcu¬ 
late x^ cit each grid point, and then find the minimum 
Xmin from which to determine Ax^ at each point and es¬ 
tablish the posterior (Eq.[^. We use a contour defined by 
^X^ = bp fo determine the boundary of the localization 
area, such that 


-f 

^ JAx^(e,<P)<bp 


e-Ax'(®>'^)/2d0d0 = p%, (17) 


where n normalizes the probability to 1. We vary the 
Ax^ value of the boundary, bp, until the integral in Eq. 


17 


reaches the desired confidence level p%. If the posterior 
is Gaussiam we can calculate bp analytic ally . The sky 
prior (Eq. 0 and the marginalization (Eq. 161 break the 
Gaussianity. However, we still use contours of constant 
Ax^ as the boundaries of our confidence level, where we 
explicitly compute the probability incorporated within 
each contour rather than assuming the an alytic value 
appro priate for true Gaussian distributions (Press et al. 
2007|). 

In summary, for a given event we use Eqs.0[^ and[T^ 
to estimate the localization area on the sky. As a test, we 
apply our lo calization algorithm to the simulated source 
catalog from I Singer et al. (|2014|), and compare our local¬ 
ization areas to those tourid from the BAYESTAR algo¬ 
rithm. As shown in Fig. we find excellent consistency. 


^ We grid the sky uniformly in (cos 6,0) with dcos9 = dip = 
0.005 (-^0.3 degree) for the two detector case. We use a finer grid- 
ding in the 3 detectors case, dcos9 = dp = 0.0025 (^0.14 degree) 
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Fig. 4. — Integrated confidence level versu s the actual captured 
percentage. The samples were taken from [Singer et al.| l |2014| | 
and localized with LIGO Hanford+Livingston\tlL,) and tlL/4-Virgo 
(HLV). The black line is the x = y line (perfect consistency). We 
find that the fraction of binaries found within our sky localization 
contours is in very good agreement with the predicted percentage. 



Fig. 5. — Ratio between the median localization area for the 
loudest event out of N detections and that for all events (N = 1). 
The SNR of the loudest events can be estimated from Eq. |18| For 3 
detectors (HLV) the areas roughly scale as SNR“^ (index=-2.04), 
while for 2 detectors (HL) the scaling is shallower (index=-1.63). 


We also explore the contribution to the localization de¬ 
termination arising from incorporating SNR and phase 
information in addition to timing, finding that the re¬ 
sulting localization areas decrease significantly (Fig. |^. 
In order to test our Bayesian posterior, we determine 
the fraction of binaries which fall within our localization 
areas, and compare this to the predicted percentages. 
Fig. a shows that our formalism is self consistent. In 
order to further study the properties of our low-latency 
localization algorithm, we simulate a population of bi¬ 
nary neutron st ar. We follow the Mo nte Carlo approach 
we presented in Chen and Ho lz (2014), utilizing the same 
noise curve as Singer et al. (|2014). h'or each binary we 
calculate the true SINK and time-of-arrival in each de¬ 
tector. We then add in Gaussian noise for each of these 
quantities, mimicking the actual measurements (ignoring 
glitches; we have assumed a threshold SNR of pth = 12 
to mimimize the importance of these). Our resulting 
distribution of localization areas is comparable to that 
in Figs. [2 andindicating that our Monte Carlo sources 
and pipwines are consistent with that of |Singer et al.| 

( [Md] ). 


3. LOUDEST EVENT 

As shown in Fig. and there is a broad distribution 
of localization areas for any given set of sources. This 
is related to the broad distribution of measured SNRs 
of the sources, which is in turn related to the distribu¬ 
tion of distances to the sources. If not all GW triggers 
can be followed, the louder (higher SNR) events offer 
the opportunity for better localization and more efficient 
follow-up. In Fig. we show the median localization 
area of the loudest source out of N events. The N = 1 
case corresponds to the full distribution of localization 
areas for all events (also seen in Fig. [^. From Fisher 
matrix calculations the localization area is expected to 
scale as Area~ SNR“^. This is true for 3 or more detec¬ 
tors, while for 2 detectors the sky locations are generally 
poorly constrained within a broad timing ring and the 
localization area will scale more slowly with improving 
SNR. In our loudest event paper (|Chen and Holz 20141 
we show that the distribution of detected SiNRs follows 


a universal distribution. We also derive the universal 
loudest event distribution, which depends upon the de¬ 
tection threshold, pth, and the number of events, N. The 
cumulative probability is: 


Fp = 1 - 

/-'max I 


N 


Pth 

Pmax 


(18) 


In (Chen and Holz 2014) we pointed out that Eq. 18 is 
independent of the detector network, d etec tor sensitivity, 
and all properties of the sources. Eq. can applied to 
all situations (e.g. 2 or 3 detectors operating at differ¬ 
ent sensitivities looking for BNS, NS-BH, or supernovae), 
and we can predict the highest SNR out of N detections 
analytically. For example, the median of this distribution 
is given by: 

_ Pth 

Pmax 


(l-0.5i/^v)i/3- 


If the localization area scales as Area~ SNR“^, we would 
expect the ratio between the median area of the loudest 
out of N events, AmaxN, and the median area of all events 
[N = 1), Al, to scale as: 


A„,axN ^PmaxN^■" ^ 1 - 0.5^^^ ^ 

We fit for the area-SNR scaling index in Fig. For 
3 detectors we see that the area scales as SNR to the 
power of -2.04, while for the 2 detector case the scaling 
is shallower (-1.63) as expected. 

Eq. [T^ shows that the localization area shrinks rapidly 
as additional binaries are initially added: waiting for the 
loudest out of the first 4 events will reduce the expected 
sky area required for follow-up by a factor of 1.9 (2 detec¬ 
tors) and 2.2 (3 detectors). This improvement declines 
for larger sets of binaries; going from 4 to 10 binaries 
reduces the sky area by an additional factor of 1.6 (2 de¬ 
tectors) and 1.8 (3 detectors). This suggests that there 
may be an optimal strategy for selecting events for EM 
follow-up, based on the expected event rate and the ob¬ 
servational constraints of the follow-up facility. 
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4. GRB JOINT LOCALIZATIONS 
Short GRBs are thought to be the very bright EM 


hole coalescence ( 

Kocevski et al. 

2010 

Fong et al.||2013 

Fong and Berger 

2U13). If a si 

lort 

GRB is accompa- 


confirm the binary nature of the short GRB engine as 
well as help constrain the binary parameters. For ex¬ 
ample, the existence of gamma-ray emission would sug¬ 
gest that at least one of the members of the binary is a 
neutron star, and measurements of the inclination would 
help constrain properties of the beamed jet. In addi¬ 
tion, the shape and area of the localization of the GRB 
in gamma-rays is independent of the Advanced LIGO & 
Virgo localization, and thus combining the two localiza¬ 
tion areas can significantly constrain the sky position of 
the source. 

We would like to explore the benefit to sky localiza¬ 
tion of a prior constraint on the inclination. Let us con¬ 
sider the case where a GRB has been observed, and we 
thereby have a prior on the GRB beaming angle which 
comes into Eqs. and [T^ We randomly generate 


in Ghen and Holz ( 

2014 

). As suggested by observations, 

we assume that GJ 

TBs 

are beamed within the range 1- 

10° (ISoderberg et al. 112006 

Burrows et al. 

2006 

Coward 

et al. 2012 hong et al. 

20L 

), and we con 

hne the orien- 


within ±cosl°, ±cos5°, ±coslO° (3 samples of 1,000 
binaries each). For each beaming angle distribution, we 
localize the source by marginalizing over the correspond¬ 
ing range of inclination angles (“sint < sin0j” which 
allows for both face-on and tail-on beaming; see Table 
or over the full range of inclinations (“random i”). We 
then compare the localization area and summarize the 
results in Table Using the beaming prior only im¬ 
proves the localization by < 10%, even in the case of 
1° beaming. This is due to the unknown direction of 
rotation; a clockwise binary produces a different phase 
pattern than a counter-clockwise binary, and knowing 
that a binary is almost face-on or almost tail-on is not 
as useful as knowing that a binary is definitely one or the 
other. If we assume that we break the clockwise/counter 
clockwise degeneracy, the localization improves signifi¬ 
cantly, approaching the “known i” case. It is to be noted 
that additional measurements of phase, e.g., from Virgo, 
LIGO-India, or KAGRA, do not help to break this de¬ 
generacy. This is because a clockwise binary “above” a 
detector produces an identical waveform to a counter¬ 
clockwise binary “below”. Alternatively, if the sky posi¬ 
tion is known, the GW network can determine whether 
the binary is orbiting clockwise or counter-clockwise. 

Although a beaming prior is of limited utility in de¬ 
termining position, the contribution of an independent 
localization from EM measurements of a GRB can be sig¬ 
nificant. The Fermi Gamma-ray Burst Monitor (GBM) 
telescope has a field-of-view of 8 steradians. It observes 
^ 45 short GRBs a year, potentially serving as a trigger¬ 
ing system for binary mergers. A limitation of the Fermi 
GBM is that it can only lo calize sources with in tens to 
hundreds of square degrees ( Briggs et al.|2009 ). We con¬ 


servatively assume that Fermi GBM localizes sources to 
a 100 deg^ (Icr) 2D circular Gaussian. This independent 




known l 

sin L < sin Oj 

random i 

r 

340 deg^ 

490 deg^ 

540 deg^ 

5 ° 

340 deg2 

500 deg2 

530 deg2 

10 ° 

350 deg2 

510 deg2 

530 deg2 

10°, clockwise 


350 deg2 



TABLE 1 

HANFORDJ-LiVINGSTON median 90% LOCALIZATION AREA FOR 
BINARY NEUTRON STARS BEAMED WITHIN 9j. ThE LOCALIZATIONS 
ARE FOR THREE DIFFERENT SETS OF PRIORS: KNOWING THE EXAGT 
INCLINATION (“KNOWN t”), KNOWING THE BEAMING ANGLE 
(“sint < sin0j”) BUT NOT WHETHER IT IS GLOCKWISE OR 
COUNTER-CLOCKWISE, AND ASSUMING NO INFORMATION ABOUT 
BEAMING (“random t”). ThE LAST ROW CORRESPONDS TO 
BINARIES THAT ARE BEAMED WITHIN 10° AND REVOLVE IN THE 
CLOCKWISE DIRECTION ON THE SKY, AND WHERE THE PRIOR 
ASSUMES NOT ONLY THE BEAMING BUT ALSO THE CLOCKWISE 
DIRECTION. 

likelihood function contributes to the posterior as an ad¬ 
ditional xlermi- Fig.|^shows the significant improvement 
in source localization which results from joint detections 
from Fermi and advanced GW detectors. The median of 
the 90% confidence 2cr) area for joint detection lies 
around 120 deg^ even with only two GW detectors, as 
compared with 460 deg^ and 560 deg^ for the GRB and 
GW localizations, respectively. In this case the beaming 
priors make no difference whatsoever in the localization. 

5. DISCUSSION 

We have developed a low latency localization algorithm 
for ground-based GW detector networks. Our approach 
requires the timing, phase, and SNR of the sources (pro¬ 
vided by the detection pipeline), and within a fraction 
of a second produces a 2D localization. We have imple¬ 
mented our algorithm for the case of non-spinning com¬ 
pact binary coalescences, assuming representative output 
from the detection pipeline (GSTLAL); this can in prin¬ 
ciple be generalized to other (preferably modeled) burst 
sources. Although the sky location is degenerate with 
other physical properties (e.g. masses, distances, and 
spins of the binaries), our algorithm avoids the need to 
estimate these paramters by instead utilizing the differ¬ 
ence in arrival times, difference in phases, and the ratio of 
measured SNRs. This increases the simplicity and speed 
of localization, while reducing potential errors. Imprecise 
waveforms can result in errors in the SNR measurements, 
but the ratio of SNRs between detectors is less sensitive 
to these errors. By taking the SNR ratio, the parame¬ 
ter space is highly reduced. For binary mergers, the two 
relevant unknowns are the binary orientation and incli¬ 
nation; it is computationally inexpensive to pre-calculate 
the marginalization of these over the entire sky The 
timing and sky priors can also be pre-gridded and pre¬ 
calculated, simplifying the localization algorithm to a 
determination using a look-up table combined with direct 
summation. Once the tables are loaded into memory, the 
entire algorithm takes < 1 second on a laptop (using a 
Python 2 script running on one thread on a 2.3 Ghz Intel 

We grid the binary inclination and orientation into 50 points 
each (2,500 points in total), and calculate the marginalization at 
each grid point. We assume that the noise curve is known and 
fixed, allowing pre-calculation of the I 7 term. Changes in the noise 
curves can be captured by the detection pipeline, and incorporated 
into our algorithm. The ratio of I 7 terms, which is the relevant 
quantity since we are concerned with the ratio of SNR, is insensitive 
to changing total mass for M < 10 Mq. 
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Fig. 6. — Cumulative distribution of sky localization for joint 
GW and Fermi GBM detections of short GRBs. This shows sig¬ 
nificant improvement, as compared with GW detection alone (Ta- 
ble0. 


Core i7). It would be straightforward to parallelize the 
algorithm (e.g., by having each core analyze a fraction 
of the sky), leading to sky localization within 0.1 sec¬ 
onds or better. Th ere is ongoing work to speed up the 
detection pipeline (Q. Chu 2015), including the possi¬ 
bility of event detection before merger. This capability 
is only of use when coupled with rapid localization; our 
prompt and accurate algorithm offers the possibility of 
real-time localization contemporaneous with the inspiral 
and merger of the binary. It is to be noted that our algo¬ 
rithm can be generalized to calculate distance as well as 
sky position, allowing for the possibility of low-latency 
3D localization. 

In Fig. we show the cumulative distribution of 50% 
localization areas for the first two years of advanced 
LIGO and Virgo network operation. Although the lo¬ 
calization is significantly improved with the addition of 
Virgo (c.f. Fig. [^, these represent a minority 25%), 
since most events detected by this network have a Virgo 
SNR below the detection threshold (pv < 4) due to 
Virgo’s l ower sensitivity and mi s-aligned antenna power 
pattern (Chen and Holz 2014). Unless the detection 
pipeline produces intorniation for these sub-threshold 
events, the low-latency localization will be based on in¬ 
formation from only two detectors (i.e., H and L), fol¬ 
lowed by better localization resulting from the full pa¬ 
rameter est imation pipeline an alyzing the sub-threshold 


waveforms ([Singer et al. 


2014); we note that there is a 
potential remedy to this discussed in Section IX of Singer 
and Price (2015). For sources detected with only 2 de¬ 
tectors, the localization of events with the 2016 network 


is potential ly worse than those detected in 2015 (Singer 


et al. 


2014). This is because the effective bandwidth ot 


the 2016 noise curve may be narrower than the 2015 case 
leading to inferior localization for sources with the same 
measured SNRs (e.g., SNR=12). Of course, the 2016 
network would detect sources from the 2015 network at 
higher SNR, and therefore would do a superior job of 
localizing identical sources. Occasionally the event will 
trigger only H-l-V or L-l-V. In these rare situations the 
localization will generally be worse than H-l-L, not only 
because Virgo operates at lower sensitivity, but also be¬ 
cause the Virgo SNR distribution skews to lower values 


when coupled with an H or L detection (Chen and Holz 


2014). 

We explore the dependence of localization on SNR, 
confirming that high SNR events lead to better local¬ 
ization. Since there is a universal distribution of SNR, 
we can use this to estimate the benefit of waiting for loud 
events. Fig. shows that waiting for the loudest out of 
the first 4 events shrinks the median localization area by 
a factor of almost 2. 

Fig. in shows that combining an EM prior on the sky 
location of a source with GW localization greatly im¬ 
proves the resulting localization. This is because the 
shapes of the two localization regions are ver different 
and completely uncorrelated. Fermi GBM detection of a 
short CRB would then greatly enhance the localization 
resulting from a GW network. One might have thought 
that such a detection would also place useful prior con¬ 
straints on the beaming angle (i.e., gamma-ray detection 
implies that the binary is face-on). However, Table 
and Fig. demonstrate that the beaming prior makes 
only a marginal improvement in localization because of 
the degeneracy between face-on and tail-on binaries. 

We have developed a simple, fast, and accurate algo¬ 
rithm for the localization of GW sources. This algorithm, 
when coupled with a fast detection pipeline, can play an 
important role in the development of multi-messenger as¬ 
tronomy. 
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